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Abstract 

I outline a perturbative QCD approach to the analysis of the deeply virtual Compton 
scattering process 7*p — > 7p' in the limit of vanishing momentum transfer t = (p' — p)'^- 
The DVCS amplitude in this limit exhibits a scaling behavior described by two-argument 
distributions F{x, y) which specify the fractions of the initial momentum p and the momen- 
tum transfer r = p' — p carried by the constituents of the nucleon. The kernel R{x,y; ^,r]) 
governing the evolution of the double distributions F{x, y) has a remarkable property: it 
produces the GLAPD evolution kernel P{x/$,) when integrated over y and reduces to the 
Brodsky-Lepage evolution kernel V{y, rf) after the x-integration. This property is used to 
construct the solution of the one-loop evolution equation for the flavor non-singlet part of 
the double quark distribution. 



^Also Laboratory of Theoretical Physics, JINR, Dubna, Russian Federation 



1. Introduction. Recently, X. Ji suggested to use the deeply virtual Compton scattering 
(DVCS) to get information about some parton distribution functions inaccessible in standard 
inclusive measurements. He considers the non-forward light-cone matrix elements 

(P _ rl2 I g(-An/2){l, -f^}^^q{Xn/2) \ P + r/2), (1) 

which appear in the lowest-order pQCD contribution to the DVCS amplitude (r is the momentum 
transfer and n a light-like 4- vector) and parameterizes them using the functions H{x,r'^,r ■ n), 
etc., with X being the Fourier conjugate variable to A. He observes that, in the r —>■ limit, the 
matrix element (0) defines the usual distribution functions like /(x), gi{x) and proposes that the 
DVCS process can be used to get information about such functions. Since the kinematics of the 
DVCS requires that r ^ 0, the = t limit can be accessed only by extrapolating the small-t 
data to t = 0. The limit r — looks even more tricky since (r • n) =0 formally corresponds to 
vanishing of the Bjorken variable xbj- Anyway, as emphasized by Ji [Q], the DVCS amplitude has 
a scaling behavior in the region of small t and fixed XBj which makes it a very interesting object 
on its own ground. 

In this letter, I briefly describe an alternative pQCD formalism for the analysis of the DVCS 
amplitude in the limit when t ^ and xbj is fixed. My main point is that, to construct a 
consistent pQCD picture for the scaling limit of DVCS, one should treat the initial momentum 
p and the momentum transfer r on equal footing by introducing double distributions F{x,y), 
which specify the fractions of p and r, resp., carried by the constituents of the nucleonQ. These 
distributions have hybrid properties: they look like distribution functions with respect to x and like 
distribution amplitudes with respect to y. Writing the matrix element of a composite operator 
in terms of the double distributions is the first step in developing the pQCD parton picture 
for the DVCS. The next step is to take into account the logarithmic scaling violation. To this 
end, I write down the evolution equation for the simplest case of the flavor non-singlet part of 
the double quark distribution F'^^ {x,y; fi). As one may expect from the preceding description, 
the relevant evolution kernel R{x,y] ^,1]) has a remarkable property: it produces the GLAPD 
evolution kernel P^^{x/^) 0, 0, ^ when integrated over y, while integrating R{x,y; ^,7]) over x 
gives the expression coinciding with the Brodsky-Lepage evolution kernel V{y,r]) for the pion 
distribution amplitude^. Using these properties of the kernel, I construct the solution of the one- 
loop evolution equation for the flavor non-singlet part of the double distribution. I also discuss 
the infrared sensitivity of the DVCS amplitude due to the presence of the light-like momentum q 
and its implications for the structure of the simplest radiative and higher-twist corrections to the 
leading- twist result. 

2. Double distributions. The kinematics of the amplitude of the process 7*p 'jp' can be 
most conveniently specified by the initial nucleon momentum p, the momentum transfer r = p — p' 
and the momentum q of the final photon. Since = 0, it is natural to use q as one of the basic 

found no advantage in introducing the average nucleon momentum P = (p + p')/2 like in eq. (|l]) 
■^Originally, the evolution equation for the pion distribution amplitude in QCD was derived and solved in ref.|^], 
where the anomalous dimension matrix Z„fc was used instead of V{x,y) (see also 0). 
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light-cone (Sudakov) 4-vectors. In the scahng hmit, the invariant momentum transfer t = r"^ and 
the square of the proton mass = can be neglected compared to the virtuality — = (g — r)^ 
of the initial photon and the energy invariant p-q = nipV. Thus, we set = and = 0, and use p 
as another basic light-cone 4- vector. Furthermore, in this limit, the requirement p'^ = (p+r)^ = p^ 
reduces to the condition p ■ r = which can be satisfied only if the two lightlike momenta 
p and r are proportional to each other: r = C^p, where C, coincides with the Bjorken variable 
C = = Q^f^ip ■ <?) which satisfies the constraint < xbj < 1- 




P p-r P P-f 

a) b) 

Figure 1: Handbag diagrams contributing into the DVCS amplitude. The lower blob corresponds 
to double quark distributions F{x,y), G{x,y). 

The leading contribution in the large-Q^, fixed-x Bj, t = limit is given by the handbag 
diagrams shown on Fig.|l|, in which the long-distance dynamics is described by matrix elements 
like 

{p-r\MOht^E{0,z;A)Mz)\p) and {p - r IMOHlf.EiO, z; A)^^) \p), (2) 

where, for Fig.^, z is the coordinate of the virtual photon vertex and E{0, z; A) is the usual 
P-exponential of the gluonic A-field along the straight line connecting and z. 

Though the momenta p and r are proportional to each other r = (p, to construct an adequate 
QCD parton picture, one should make a clear distinction between them. The basic reason is that 
p and r specify the momentum flow in two different channels. For r = 0, the net momentum 
flows only in the s-channel and the total momentum entering into the composite operator vertex 
is zero. In this case, the matrix element coincides with the standard distribution function. The 
partons entering the composite vertex then carry the fractions of the initial proton momentum 
(— 1 < Xi < 1). When x is negative, the parton is interpreted as belonging to the final state and 
Xi is redefined to secure that the integral always runs over the segment < x < 1. In this parton 
picture, the spectators take the remaining momentum (1 — x)p. On the other hand, if the total 
momentum flowing through the composite vertex is r, the matrix element has the structure of the 
distribution amplitude in which the momentum r splits into the fractions yr and {1 — y)r = yr 
carried by the quark fields attached to that vertex. In a combined situation, when both p and r are 
nonzero, the initial quark has the momentum xp + yr, while the final one carries the momentum 
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xp — yr. In more formal terms, this corresponds to the following parameterization of the light-cone 
matrix elements 

{p-r\MO)zE{0,z;A)M^)\p)U2^o = u{p-r)zu{p) /' /' (e-'^^P'^-^y^'''^ F,{x,y) 

Jo Jo ^ 

_ e«(p-)-^5(-)F^(a;, y)) e{x + y< l)dy dx, (3) 

{p-r\ M0)j5SEi0, z- A),Pa{z) \ p) |,2=o = u{p - r)-i,zu{p) (e-'<P^^-'y^'''^Ga{x, y) 

+ e'''^P'^-'y^''^Ga{x, y)) e{x + y< l)dy dx, (4) 

where z = j^z^ and u{p — r),u{p) are the Dirac spinors for the nucleon. 

Though we arrived at the matrix elements (^), (^) in the context of the scaling limit of the 
DVCS amplitude, they accumulate a process-independent information and, hence, have a quite 
general nature. The coefficient of proportionality between p and r is then just the parameter 
characterizing the "asymmetry" of the matrix elements. The fact that, in our case, ( coincides 
with the Bjorken variable is specific for the DVCS amplitude. The most non-trivial feature implied 
by the representation (^) is the absence of the ^-dependence in the double distributions Fa{x,y) 
and Ga{x,y). Using the methods developed in ref.0, one can prove that this property and the 
spectral constraints x>0,y>0,x + y<l hold for any Feynman diagram. As a result, both the 
initial active quark and the spectators carry positive fractions of the light-cone "plus" momentum 
p: x + (y for the active quark and x — (y= {l — x — y) + {l — ()y for the spectators. Note, however, 
that the fraction of the initial momentum p carried by the quark going out of the composite vertex 
is given by x — y( and it may take both positive and negative values. At first sight, this result 
contradicts the intuition based on the infinite momentum frame picture. Recall, however, that all 
the fractions are positive only when the "plus" component of the momentum transfer r vanishes, 
which is not the case here, since r has only the "plus" component. 

Taking the limit r = gives the matrix element defining the parton distribution functions 
fa{x), fa{x) (or ga{x), Qa^x)) . This obscrvatiou results in the following reduction formulas for the 
double distributions F{x,y),G{x,y): 

/ Fa{x,y)dy = faix) , / Ga{x,y) dy = ga{x). (5) 
JO Jo 

3. Leading-order contribution. Using the parameterization for the matrix elements given above, 
we get a parton-type representation for the handbag contributions to the DVCS amplitude: 

T^^'ip^q^r) = - ^{p^q'' +p^q^)\ J] (T^(C) + T^(C)) 
\ V ' q J a 

Me^-'p^q, ^^^^ E el (T^(C) - Tj(C)) , (6) 
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where q = 7^?^^, the factor a/1 — C, comes from u{p') = y/1 — Cu{p) and Ty{(), T^(^) are the 
invariant amphtudes depending on the scahng variable (: 

The terms containing l/{x — (y + ie) generate the imaginary part: 
1 , /•! /-i 



-ImT^(C)=/ / 5(a;-Cy)Fa(a:,y)^^(a; + y < l)rfxrf|/ 

TT JO Jo 

f Fa{x,l-x/Odx= f' Fa{yC,y)dy, (9) 
JO Jo 



-I' 



with a similar expression for ImT|((^). Because of the integration remaining in eq.(P), the relation 
between lm.T{Q and the double distributions Fa{x,y) is not as direct as in the case of forward 
virtual Compton amplitude, the imaginary part of which is just given by distribution functions 
faiO- Note, that the ?/-integral in eq.(|^) is different from that in the reduction formula (H), i.e., 
though 

$a(C)= f'FM,y)dy (10) 

Jo 

is a function of the Bjorken variable (, it does not coincide with /a(C)- 

To get the real part of the l/{x — (y + ie) terms, one should use the principal value prescription, 
i.e., ReT(C) is related to Fa{x,y) through two integrations. 

4- Evolution equation. In QCD, the limit —>■ for the matrix elements in eq.(^ is singu- 
lar. As a result, in perturbation theory, the distribution F{x, y) contains logarithmic ultraviolet 
divergences which can be removed in a standard way by applying the i?-operation characterized 
by some subtraction scale /i: F{x, y) F{x, y; /i). The /i-dependence of F{x, y; jj) is governed by 
the evolution equation 

( d d\ 

\dfl ^ ^^^"^dg ) ^' ^ Jo Jo ^' ^' ^' ^' ^'^'^^ 

Since the integration over y converts F{x,y) into the parton distribution function f{x), whose 
evolution is governed by the GLAP equation 

+ f^(9)-^^ fix; /i) = £ jP{x/^; g)f{i; fi)d^, (12) 

the kernel R{x, y; ^, rj; g) must have the property 

/'""^ R{x, y; i, v; g)dy = -P{xli). (13) 
Jo 4 
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For a similar reason, integrating R{x, y; ^, rj; g) over x one should get the Brodsky-Lepage kernel: 



R{x,y;^,V;9)dx = V{y,ri;g). (14) 
Explicit one-loop calculations give the following result for the gg-component of the kernel: 

y; g) = '^Cf-^ |^(0 < x/i < min{|//r7, y/n]) ^imv - v) (15) 

6(0 < x/i < l)x/i 
(1 - ^/O 



+ 



Jo I — z J 



Here the last (formally divergent) term, as usual, provides the regularization for the singularities 
of the kernel for x = C, (or y = rj). Note that, in the Feynman gauge, the first line of eq. (|T5|) 
corresponds to operators with ordinary derivatives d'^ while the second one results from the d'^ — >■ 
jju _ Qv _ ig/^^ change. It is easy to verify that the kernel Rqq{x, y; ^, rj; g) has the property that 
X + y < 1 ii ^ + rj < 1. Using our expression for Rgq{x, y; ^, rj; g) and exphcit forms of the Pgg{x/^) 
and V{y,ri) kernels 



TT 



l + z' 
1 - z 



v{y.v) = ^cMy- 



1 + 



1 



n-y 



e{y < T]) 



1 + 



1 



e{y > rj) 



(16) 
(17) 



y-v_ 

one can check that Rqq{x,y; ^,1]; g) 



(where "+" denotes the standard "plus" regularization | 
satisfies the reduction formulas (|T3D and (|14|). 

Note also that the flavor non-singlet light-cone operator 0^^{z, 0) = '^l!{z)X°'z^lJ{0) satisfies the 
Balitsky-Braun evolution equation H| 



du I K{u,v)0^^{uz,vz)dv 



This means that our kernel Rgq{x, y; ^, rj; g) should be related to the i^-kernel by 

Rgq{x,y;i,r];g) = ^K{y + r]x/^,y - f]x/^). 
Indeed, taking the explicit form of K{u, v) from ref. 00 



K{u, v) 



TC 



Cf[1 + 6{u)[v/v]+ + 6{v)[u/u]+ - -6{u)6{v] 



(19) 



(20) 



and combining it with eq. ([T9|) , one immediately obtains our expression ( |T5|) for Rgg{x,y;^,ri;g). 

To solve the evolution equation, I propose to combine the standard methods used to find 
solutions of the underlying GLAP and Brodsky-Lepage evolution equations. Here, we will consider 

^In the definition adopted in ref. K{u, v) has the opposite sign. 
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only the simplest case, i.e., the evolution equation for the flavor-nonsinglet component in which 
qg, gq and gg kernels do not contribute. Recall first that, to solve the GLAP equation, one should 
consider the moments with respect to x. Integrating x'^Rgq{x,y; ^,1]; g) over x and utilizing the 
property Rqq{x,y;^,r]; g) = Rqq{x/^,y; 1,7]; g)/^, we get 



(4 ^'"''S 

where Fn{y; jj) is the nth x-moment of F{x, y; /i) 



Fn{y]^) 



x'^F{x, y; fj,)dx 



(21) 



(22) 



and the kernel Rn{y, d) is given by 

n+l 



a. 



Rn{y,v;9) = —Cf 

TT 



y 



1 1 

+ 



n + 1 f] " y 



e{y < v) 



' \ n+l 



1 1 

+ 



n + l y ~ "fl 



e{y > r,) 



1 — z 



dz\. (23) 



It is straightforward to check that Rn{y,rj;g) has the property 

Rnijj, m g)wn{ri) = Rnij], y; 9)wn{,y), 

where Wniy) = {yyY^^. Hence, the eigenfunctions of Rn{y,ri;g) are orthogonal with the weight 
Wn{y) = {yyY^^ 1 i-^., they are proportional to the Gegenbauer polynomials C^^^^'^{y — y) (cf.[^ 
|I0]). Now, we can write the general solution of the evolution equation 



Fn{y: fi) = {yyT^' E AnkCT'"\y - y) [log(/x/A)] 

fc=0 



(24) 



where /?o = 11 ~|^/ is the lowest coefficient of the QCD /3-function and the anomalous dimensions 
7^"'' are related to the eigenvalues of the kernel Rniy-, V'l 9)'- 



(n) 

Tk 



1 -r fc+n+1 1 

+ 2 V - 

2 (n + A; + l)(n + A; + 2) fr'^ J 



(25) 



They coincide with the standard non-singlet anomalous dimensions IjTT], 0: 7^"^ = 'jn+k+i- 
Note, that 7q'^^ = 0, while all other anomalous dimensions are positive. Hence, in the formal 
/i — s> 00 limit, we have Fo{y, /i — > 00) ~ yy and Fn{y, fi 00) = for all n > 1. This means that 

F{x,y;n-^ 00) ~ 6{x)yy, 

i.e., in each of its variables, the limiting function F{x,y;fi 00) acquires the characteristic 
asymptotic form dictated by the nature of the variable: 6{x) is specific for the distribution func- 
tions [jiT], |12[, while the yy-foim is the asymptotic shape for the lowest-twist two-body distribution 
amplitudes M, M. 
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5. Infrared sensitivity and hard gluon exchange corrections. The structure of the leading term 
of the DVCS amphtude T(C, Q"^) is very similar to that of the forward virtual Compton scattering 
amplitude Tf{uj, Q"^) which is the starting point of the DIS analysis. The major difference between 
the two amplitudes is that one of the photons in the DVCS amplitude is real. As a result, in 
higher orders, the coefficient functions C{(, Q^, q^) of the formal operator product expansion of two 
electromagnetic currents may be singular (non-analytic) in the limit — > 0. These singularities 
are related to the possibility of a long-distance propagation in the g-channel. In fact, the long- 
distance sensitivity simply means that the relevant contribution is non-calculable in perturbation 
theory, and one should describe/parameterize it by introducing the distribution amplitude for the 
real photon. 




P p-r P P-r 

a) b) 

Figure 2: Simplest hard gluon exchange contributions to the DVCS amplitude. The upper 
blob corresponds to the photon distribution amplitude (Pj{u) and the lower one to double quark 
distributions F{x, y), G{x, y). 

Two simplest contributions of this type are shown in fig.^j. They are analogous to the hard 
gluon exchange diagrams for the pion form factor. The basic common feature is that the large 
virtuality flow bypasses the real photon vertex, and the relevant contribution factorizes into the 
product of the hard scattering (short-distance) amplitude formed by quark and gluon propagators 
and two long-distance parts. In our case, one long-distance part is given by the double distribution 
F{x,y) (or G{x,y)) and the other by the photon distribution amplitude (p-y{u). The latter can be 
understood as the probability amplitude to obtain the photon with momentum q from a collinear 
qq pair, with quarks carrying the momenta uq and (1 — u)q. The w^-moments of the lowest-twist 
function ^p^{u) are related to the difference between the "exact" and perturbative versions of the 
correlator 

^,uu,...uM = / e^''{Q\T{J^{z) 4(0)7.^., ■ ■■D,^i,Mmd'z (26) 

of the electromagnetic current J^{z) with a composite operator containing n covariant derivatives. 
The perturbative amplitude for the qq ^ 'y transition is taken into account in the pQCD radiative 
correction to the coefficient function (c/. [0). By a simple counting of the propagators, one can 
easily find out that the contribution of Fig.^ behaves like as/ in the large-Q^ limit. Furthermore, 
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due to the EM current conservation, this correlator is proportional to the tensor structure Qfj^u^f ~ 
g^g^. Since, for a real photon, = and (ge) = (e^ being the photon's polarization), the relevant 
contribution vanishes. Hence, the hard gluon exchange corrections to the DVCS amplitude are 
rather strongly suppressed. 

At the leading-power level 0{Q^), the absence of the non-analytic ln(g^)-terms in the one- loop 
coefficient function contribution for the A-part of the amplitude can be read off from existing 



results for the correction to the 77* vr amplitude [|14|, |T5|, jl 



6. Conclusions. In this paper, I formulated the basics of the pQCD approach for studying 
the deeply virtual Compton scattering amplitude in the scaling limit. The essential point is 
that information about the long-distance dynamics in this case is accumulated in the double 
distributions y). I studied the evolution properties of y) and also showed that ?/)'s 
are related to the standard distribution functions f{x) through the reduction formulas (^. An 
interesting problem for the future studies is to understand a possible interplay between the x- and 
^-dependencies of the double distributions and construct phenomeno logical models for F{x, y) . 
The simplest idea is to try a factorized ansatz F{x, y) = f{x)g{y). However, the explicit expression 
for the evolution kernel R{x, y; ^, t]) (which satisfies a similar reduction formula (|13|) ) suggests 
that the structure of F{x,y) is not necessarily as simple as that. 

In the experimental aspect, as emphasized by Ji [0, a continuous electron beam accelerator 
with an energy 15 - 30 GeV (like proposed ELFE) may be an ideal place to study the scaling 
limit of the DVCS. It is also worth studying the possibility of observing the first signatures of the 
scaling behavior of DVCS at CEBAF, especially at upgraded energies. 

In a forthcoming paper [|1^, I discuss the gluonic double distributions Fg{x,y) which play a 
crucial role in the perturbative QCD approach to hard diffractive electroproduction processes like 
7*p — > p'p. 
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